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Abstract 

An exact two-dimensional rotation-strain model describing the motion of 
Hookean incompressible viscoelastic materials is constructed by the polar decom- 
position of the deformation tensor. The global existence of classical solutions is 
proved under the smallness assumptions only on the size of initial strain tensor. 
The proof of global existence utilizes the weak dissipative mechanism of motion, 
which is revealed by passing the partial dissipation to the whole system. 

Keyword: Rotation-strain model, viscoelastic fluids, partial dissipation, complex flu- 
ids. 

1 Introduction 

Viscoelastic materials include a wide range of fluids with elastic properties, as well 
as solids with fluid properties. Let us start with the standard description of general 
mechanical evolutions to introduce some notations and definitions. Any deformation 
can be represented by a flow map (particle trajectory) x(t, X), < t < T, a time- 
dependent family of orientation-preserving diffeomorphisms. X is the original labeling 
(Lagrangian coordinate) of the particles. X is also referred to as the material/reference 
coordinate, x is then the observer's (Eulerian) coordinate. In general, the velocity field 
u(t, x) is defined as the time derivative of the flow map, which is a vector field defined 
on the Eulerian coordinate: 

u(t, x) = x t (t, X(t, x)) . 
The deformation tensor F in Lagrangian coordinate is defined by 

F 3 (t,X) = ^-(t,X). 

When one uses the Eulerian description, the deformation gradient tensor F in spatial 
coordinate is defined by 

F(t,x) = F(t,X(t,x)). 
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Recently, an extensively studied system modelling the motion of isotropic incom- 
pressible viscoelastic materials of Oldroyd-B type takes: 



V -u= 0, 

u t + u- Vu + Vp = fiAu + V-(^PF T ), (i.i) 
F t + u ■ VF = VuF, 

where p(t, x) is the pressure, \i is the viscosity and W(F) is the elastic energy function. 
The global existence of classical solutions near equilibrium to (11. ip in the 2D Hookean 
case was established by Lin, Liu and Zhang [17] (see also Lei and Zhou [15] via the 
incompressible limit method). Later on, Lei, Liu and Zhou [H] proved the global 
existence of classical small solutions to system ( 11. ip in both 2D and 3D cases (see also 
Chen and Zhang [3J). Very recently, an improved result was obtained by Kessenich [9] 
by removing the dependence of the smallness of the initial data on the viscosity via 
the hyperbolic energy method. 

To study the different contributions of strain and rotation parts, Friedrich [5] ob- 
served that the smallness of the strain in nonlinear elasticity can be realized through 
the polar decomposition of the deformation tensor. John P, [7J showed that no point- 
wise estimate for rotations in terms of strains can exist even in the case of small strain. 
In the work of Friedrich and John, no PDE is involved. From the PDE point of view, 
Liu and Walkington [19J considered approximating systems resulting from the special 
linearization of the original system with respect to the strain. Lei, Liu and Zhou [13] 
constructed a 2D rotation-strain viscoelastic model and proved the global existence 
their classical solutions with small strain. One shortcoming of the model in [T3] is 
that the dynamics of strain and rotation are not equivalent to that of the deformation 
tensor. 

One main concern of this paper is to construct an exact rotation-strain model for 
motions of Hookean viscoelastic materials of Oldroyd-B type, which takes 

V-m = 0, 

u t + u ■ Vu + Vp = /iAu + V • {VV T ) + 2V • V, 

V t + u-VV = D(u) + l(VuV + VVu T ) 

+ \lo X2 {u){VA - AV) - \>y(VA - AV), (1 ' 2) 



vM + ,w^>-p: v v ;) + vv 

Here V is the strain part of the deformation tensor, 9 is the rotation angle and 7 is a 
nonlinear term. The derivation of this model will be presented in section 2. It seems 
that the velocity field u is the only dissipative unknown. We will prove the global 
existence of classical solutions to ( j 1.2ft for small initial strain V and initial gradient of 
9 (9 may be arbitrarily large). The proof utilizes the weak dissipative mechanism of 
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system ( 11.21) by passing the dissipation of the velocity field u to the strain tensor V and 
the gradient of the angle variable V0. Note that the equations for u and V have already 
formed an closed dynamical system. However, our proof needs the use of the equation 
for 9 to carry out the weak dissipative mechanism of the whole system. Without the 
help of the equation for 9 and the underlying physical relationship between V and 9 
(see Lemma 13731) . the corresponding analysis issue is still open. 

At last, let us also cite some related results on elastodynamics when the viscosity 
fi is zero. The global existence of 3D compressible elastodynamics was established by 
Agemi pQ and Sideris [20] independently under the null condition and small initial data. 
An earlier almost global existence result can be founded in John and Klainerman [8] (see 
also Klainerman and Sideris [TT]). Hyperbolic energy method involving Klainerman's 
vector field developed by Klainerman [10] (see also Christodoulou [I]) plays an essential 
role in [HI I2Q] . The incompressible case was later on studied by Sideris and Thomases 
in J2TJ [221 123]. A natural question to ask is that what happens for the solutions of 
( 11. ip when the viscosity vanishes, which is answered very recently by Kessenich [9] in 
the case of 3D. Other results on Oldroyd-B model can be found in [21 [T2J [161 ITS] . 

The remaining part of this paper is organized as follows. In section 2 we present the 
derivation of the rotation-strain model. Then the main result of this paper is stated. 
Section 3 is devoted to exploring some of the intrinsic properties of the viscoelastic 
system. The proof of global existence theorem is completed in section 4. 



2 Derivation of the Strain-Rotation Model 

Constitutive relations generally involve the deformation tensor F. In particular, the 
polar decomposition of F is important. For any non-singular matrix F, it is easy to see 
that FF T is positive definite symmetric. Thus, there exists a unique positive definite 
symmetric matrix I + V such that (I + V) 2 = FF T , where / is the 2x2 identity matrix. 
Then the matrix R = (I + V)~ 1 F is orthogonal. Consequently the non-singular matrix 
F can be decomposed uniquely in the form 

F=(I + V)R, (2.1) 

where R is orthogonal 

RR T = I 

and I + V positive definite symmetric 

V = V T . (2.2) 

Physically, this means that the deformation is decomposed into stretching and rotation. 
Following a suggestion by K. O. Friedrichs [5], the tensor I + V is called the left stretch 
tensor, V the strain matrix and R the rotation matrix. 

Our first goal is to formulate the viscoelastic system in terms of the velocity u 
and the strain V. We first look at the momentum equation, the second equation in 
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(II. lp . We restrict our discussions below to a special elastic energy functional of the 
Hookean linear elasticity W(F) = ^\F\ 2 , it does not reduce the essential difficulties for 
mathematical analysis. Indeed, all the results we describe here can easily be generalized 
to a more general class of elastic energy functionals of the deformation tensor F (see 
|14j). In terms of the strain matrix V, the momentum equation takes the form of 

ut + u ■ Vu + Vp = fiAu + V • {VV T ) + 2V • V (2.3) 

where we used the equation ( 12. ip and ( 12. 2ft . We point out here that through this paper, 
V with no indices stands for the derivative with respect to the spatial coordinate x. 
For a matrix-valued function A, (V--B)i means VjBij, where summation over repeated 
indices will always be well understood. 

In order to get a dynamical equation for V, we plug (12.11) into the third equation 
of system ( 11. ip . Denote the material derivatives of a quantity Q by 

Q = Q t + u-VQ. 

We have 

VR+(I + V)R = Vu(I + V)R. 
Noting that R is orthogonal, we deduce that 

V + (I + V)RR T = Vu(I + V). (2.4) 
By transposing the above equality (12. 4p . it is rather easy to see 

V + RR T (I + V) = (I + V)Vu T . (2.5) 
Now let us write the orthogonal matrix R as 



R 

It is immediately verified that 
RR T = 



cos 9 — sin 6 
sin 9 cos 9 



Define 



we have 



-sin# — cos^A / cos# sin#\ a_ fO — 1 
cos 9 — sin 9 ) \ — sin 9 cos 91 \ 1 



A = (° - 1 ) , (2.6) 



RR T = A9, VA + AV = AtiV A T = A' 1 = -A. (2.7) 
Thus, it follows from (E3J), (ESJ and (EZD that 

(J + - A T (J + V)9 = Vu(I + V)-{I + V)Vu T 
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which is equivalent to 



91 = -^-±- v A[Vu(I + V)-(I + V)Vu T ] 

[2u} 12 (u) + V k uiV k2 - V fc w 2 Vfei]/ 



2 + trV 

-LU 12 (u)I + 7/, 



where 



u(u) = ]-(Vu-Vu T ) (2.8) 



is the virticity tensor and 



7 = o \ ,J trVu 12 (u) - {V k u x V k2 - V k u 2 V kl )}. (2.9) 
2 + trV 

Consequently, we obtain that 

= -u; 12 (u) + 7- (2-10) 

Next, add up ( 12. 4ft and ( 12. 5ft and then insert ( 12. 7ft and 1 12. 10p into the resulting 
equation, we obtian 

1 



V = ^[Vu(I + V) + (I + V)Vu T ] (2.11) 
1 



- |[(/ + ^)i?i? T + i2i? T (/ + V)] 

= ^ ( Vw + Vw T ) + - ( VuV + V Vu T ) 

+ ±{(I + V)A + A t (I + V)]lj 12 (u) 

- 1 - 1 [(I + V)A + A T (I + V)] 

= l - ( Vu + Vw T ) + i ( VmF + V Vu T ) 

+ ~u; 12 («)(WL - AV) - ly{VA - AV), 

where in the last equality we used (12.71) . If we denote the symmetric part of the gradient 
of velocity by 

D(u) = ^{Vu + Vu T ), (2.12) 
we can rewrite the above equation (12. lip as 

V t + u-VV = D(u) + -(VuV + VVu T ) (2.13) 
+ ~ Ul2 (u)(VA-AV)-^ 1 (VA-AV). 



In order to formulate an incompressible dynamical problem, the equations of mass 
conservation V ■ u = and balance of momentum (12.31) supplemented by the transport 
equation (12.131) have in fact made up of a closure system. However, to prove the global 
existence, we must include an equation for V# to carry out some intrinsic properties 
of the rotation-strain viscoelastic system, which will be seen more clear in section 4. 

Let 



V- = ^ v ;j (2.14) 
and then apply V ± to (12. 101) to yield 

+ u • VV^ = '-An - V V ;) + VH (2.15) 

This equation will be used to explore the weak dissipation of V#, which eventually 
shows that V is also weakly dissipative. The equations (12. 3p . (12.131) and (12.151) form 
the rotation-strain model (11.21) for viscoelastic fluids. For future use, we also need an 
equation of V8. For this purpose, we apply the gradient operator V to (12.101) to get 

ve t + u ■ we = - (l 1 ^ ' ^ - vwia(«) + v 7 . (2.16) 

We will consider the Cauchy problem or the periodic boundary- value problem of 
system ( 11.21) . The initial data takes of the form 

u(0,x) =u (x), V(0,x) = V (x), V0(O,x) = V0 o (x), xGO. (2.17) 

where Q C R 2 is two dimensional torus or the entire space. The above initial data will 
be imposed on the following constraints 

V-« = 0, 

det(/ + V ) = l, (2-18) 
V-V = A(I + V )V9 . 

The first two are just the consequences of the incompressibility and the last one is 
understood as the consistency condition for changing of variables (see Lemma 13.11 
Lemma 13.31 Remark 13.21 and Remark 13.41 in the next section). 

Remark 2.1. In our earlier work [13], the rotation matrix R satisfies a spontaneous 
but specified transport equation: 

Rt + u- VR = -ARu 12 (u) = uj{u)R. 

Consequently, it follows that 

6 t + u ■ V9 = -wia. 
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This combining the third equation of system f 1 1.1 1) gives a non-symmetric transport 
equation for V. While in the model (jl.2p . noting (12. 7p and with the aid of (I2.10p . it is 
rather easy to find that 

R t + u-VR = ARB = -ARuj 12 (u) + 7AR. (2.19) 

The last term of the right side of the above equation or the one 7 in (I2.10p repre- 
sents corrections of our previous model, which leads to that the strain matrix must 
be symmetric. However, the corrections reflect the intrinsic properties of motions of 
viscoelastic fluids and the underlying physical origins. 

For the rotation-strain model (11.21) . we will prove the following theorem: 

Theorem 2.2. Consider the Cauchy problem or the periodic initial-boundary value 
problem for the rotation- strain viscoelastic model ( 11. 21) and ( 12. 17ft with the intrinsic 
physical constraints ( 12.181) on the initial data. Then there exists a unique global classical 
solution (u, V, V8) which satisfies 



\Hh(si) + \\V\\h(si) + W Au + -V^f + n 

t 1 

+ ||V(Au + -V^)f + ||V(Au + -V • V)\\' 
if the initial data u , V G H 2 (Q) and V6* G H 1 ^) and 



|Vw" 2 



H 2 



dt < 



^ 4 (1 + /) 
C 2 (l + /i 10 )' 



holll^n) + \\ v o\\ 2 Hi(n) + \\^ 9 o\\m(n) < ^oy 

where Q C R 2 is two dimensional torus or the entire space, C and M (M > 2C 3 ) are 
big enough constants independent oft and \i. 

The proof of the above theorem relies on a local existence theorem and a priori 
energy estimates. The proof of local existence is standard and can be similarly done 
as in [T^l IT6] . thus is omitted here. Below we will only present the a priori estimate 
in Theorem r 



Remark 2.3. In the case of Q = T 2 being a two-dimensional torus, we in fact recovered 
the results in our earlier work [14|. To see this, we set 

If 

6*o = 77^77 / 9 dx. 



IT 2 



T 2 



By Poincare's inequality, we have 

||0o -0o\\l°° < C\\e -9 \\ H 2 < C\\V9 \\ H i. (2.20) 
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Thus, if we define 

= /cos^o -sin0 o \ 

ysin 9 cos 6> / 

then 

F(0,x)= [(I + V)R] t=Q = R + E 

where 

E = [R-R + VR]\ t=0 

is a small disturbance near a constant equilibrium R by (12.201) . since HWollff 1 is 
assumed to be small in Theorem 12.21 In [TJ] , we studied the general viscoelastic model 
including an equation for the deformation tensor F and proved the global existence of 
classical solutions near constant equilibrium in both two and three-dimensional spaces. 
However, our results and proofs here are still of importance even in the periodic case 
because they provide us a better understanding of the physical background of the 
system. 

3 Special Structures of the System 

In this section, we will explore some of the intrinsic properties of the viscoelastic system 
mentioned in section 2. These properties reflect the underlying physical origins of the 
problem and in the meantime, and are also essential to the proof of the global existence 
result here. 

The following lemma will be used to explore the weak dissipation of V. 

Lemma 3.1. Assume that the second equality of (12.181) is satisfied and (it, V) is the 
solution of system (11.21) . Then the following is always true: 

det(/ + ]/) = l (3.1) 

for all the latter time t > 0. In other words, we have 

trV = - det V. (3.2) 

Proof. Using the identity that for any non-singular matrix F, 

ddetF 



dF 

we have 



det FF~ 



(det(I + V)) t + u-V(det(I + V)) (3.3) 
= det{I + V){I + V)^V tJ 



s 



Return to (12. lip , we in fact have that 

V = D(u) + -{VuV + VVu T ) - ^(VA - AV)6. (3.4) 

By 1ET5J) . (E3D, we arrive at 

(det{I + V)) t + u-V(det(I + V)) (3.5) 

= trD(u) - ^tr(VA - AV)6 + trl^VuV + D ^ 

We will compute the right side of the above equality term by term. First of all, 
with the use of the first equation of (11.21) . it is rather easy to see that 

trD{u) = 0. (3.6) 

Making use of (I2.6p . we deduce that 

According to the expression of the Cauchy strain tensor (I2.12p . the third term of (13. 5ft 
can be calculated as 

tr { VuV+ {-v 12 ^) D ^} < 3 ' 8 > 

= VlUlKi + V 2 UiV l2 + V1U2V12 + V 2 -U2^22 

+ V22V1M1 - Vi 2 (Viu 2 + V 2 «i) + V u V 2 u 2 
= trVW ■ u = 0. 

Similarly, by a straightforward computation, we find that 

JPk. "^Wl (3.9) 



V22 - 


-v 12 \ , 


-V12 




/ ^22 


-v l2 


\-Vv2 





Vi 2 \ (ViUiV n + V 2 u x Vi 2 V1U1V12 + V2M1V22 

WlU 2 Vu + V 2 U 2 Vi 2 ViU 2 Vi 2 + V 2 M2^22 



V 22 (V l u 1 V 11 + VaUi^is) + Ki(Viit 2 ^i2 + V 2 u 2 V 22 ) 

- V 12 [(V lU2 Vn + V2W2V12) + (Vittiyia + V 2Ul V 22 )} 
V n V 22 \7 ■ u + V r i 2 (^2V 2 ui + VttViu 2 ) 

- ^i 2 [ViW 2 Fn + V2M1V22 + V 12 V ■ u] 
0. 
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It remains to deal with the last term of (13.51) . Recall (12.61) once more, it follows from 
a simple calculation that 




V 22 —V\2 

-Via V u 



(VA-AV)J (3.10) 

V 22 -V 12 \ ( 2V 12 V 22 - V n 
-V 12 V u J \V 22 - V n -2V 12 

= 2V 22 V 12 - 2V 12 {V 22 - V n ) - 2VuV X2 = 0. 

Finally, combining the above equalities (13.51) with (I3.6l) - (l3.10p . we arrive at 

(det(J + V)) t + u ■ V(det(/ + V)) = 0. (3.11) 

This completes the proof Lemma 13.11 

Remark 3.2. The incompressibility can be exactly represented as 

detF = l. (3.12) 

The usual incompressible condition V • u — 0, the first equation in (jl.2p . is the direct 
consequence of this identity. On the other hand, det(J + V) = 1 is also a direct 
conclusion of ( 13.121) since R is an orthogonal matrix. However, the above Lemma 
illustrates the incompressible consistency of the the system ( II. 2(1 . 

To show the global existence of the classical small strain solutions to system ( 11.21) . 
we still need to show that the pressure is a high order term, which will be seen more 
clear in the next section. To do so, we will use the following lemma to find out the 
weak dissipation of V0, which will play an important role for our proof below. 

Lemma 3.3. Assume that the third equality of f 12. 18(1 is satisfied. Then any solution 
(u,V) of the system (I1.2I) - (I2.17I) will satisfy the following identity 

V -V = A(I + V)V6 (3.13) 

for all the latter time t > 0. 

Proof. We rewrite the third equation of (11.21) as 

V t + u ■ W = D(u) + ^{VuV + VVu T ) - -{VA - AV)9, (3.14) 

and then apply V- to the resulting equation (13.141) to yield 

V-Vt + M-VV-V (3.15) 

- -VAV9 - -v ■ (VA - AV)0 

= 5 Ao + \ Avve - (vj: : wS) + \ v ■ < v " y + KV " T » 

- WAV6 - -V • {VA + AV)9 + AV ■ V6. 
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By the first equation of (II. 2p and a straightforward computation, we obtain 



v /; V1 ^ + -v ■ (Vuv + Wu 1 ) = -v ■ (v^y - vv« J ). (3.16) 

In view of (12.101) . we can rewrite the last term of (13.151) as 

AV -Ve =[u(u) + ~/A\V -V. (3.17) 
Thus, it follows from (I3TT5D . (l3TT5j) and (13371) that 

V -V t + u-VV -V (3.18) 

= ^An + ^Avve + [u(u) + jA] v ■ y - ^vAve 

+ ■ (Vu^ - VVu T ) - • (VA + AV)6. 

On the other hand, combining (12. 161) with (13.141) . and splitting some terms of the 
resulting expression, we conclude that 

d t [A(I + V)V9] + u ■ V[A(I + V)V6] (3.19) 
= AD(u)V9 + -A(VuV + VVu T )Ve - ^(AVA + V)9V6 



1 1 ..—A i 



+ A(/ + y) V0 - A(/ + v) 



= Ave + -Awe + -Awe + v# + -A(vuv + y w )w 

In view of (12.101) and the first equation of system (jl.2p . we obtain, after a straightfor- 
ward calculation, 

AV6 = -Au + AV'y. (3.20) 



It is rather easy to see that 

— A I 

v 2 u ■ ve 



AD(u)Ve -A[Y} U ' Y_ 6 A = u{u)AVe (3.21) 



and 



r \r\-!„.T\ 



A(VuV + VVu 1 )Ve - AV I ^ L u ^ Q 1 = -A(VwF - )V0. (3.22) 



Using (|22D and ([2TUD . we find that 



--(,4FA + F)0V0 (3.23) 
= -^A(VA + AV)eve - {-to 12 {u) + 7 )\/w 
= -^A(VA + AV)(?V0 + Au;(w)VV0 - 7 VV0. 
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Combining fl3.19p - fl3.23l) . we have 

d t [A(I + V)V9] + u ■ V[A(I + V)V9] (3.24) 
= ^Au + ^AVV9 + uj(u)A(I + V)V9 - ^A(VA + AV)9V9 

+ AVj + -A(VuV - VVu T )V8 + -AVV6 - 7 W0 

Now, we subtract (I3.24p from (13.181) . and rearrange the right side terms of resulting 
equation. The outcome of this straightforward calculations is 

d t [V ■ V - A(I + V)V8] + u ■ V[V • V - A(I + V)V9] (3.25) 
= u{u) [V-V-A(I + V)V9] - -(VA + AV)V9 - - V • (VA + AV)9 

+ - V • (VuV - VVu T ) - AV 7 + 7AV • V 

+ [-A(VA + AV)9 - -A(VuV - VVu T ) + -yV]V9 

— u(u)[V ■ V - A(I + V) V9] - • [(VA + AV)9] 

+ ■ (VuV - VVu T ) - AV'j + 7AV ■ V 

+ [^A(VA + AV)9 - ~A(VuV - VVu T ) + jV]V9. 

Our objective is to show that the expression V- V — A(I+V) V9 satisfies some differential 
equation whose solution can be uniquely determined by the initial data. Recall (12. 7\i 
and fT2TTUj) . we get 

- -V ■ [(VA + AV)9) + ■ (VuV - VVu T ) - AVj + jAV ■ V 

= - -V • [trV A(-u 12 (u) + 7)] - AV7 + 7AV • V 

+ iv • (VuV - VVu T ) 

= -AV[trVu 12 (u)} - -AV[(trV + 2)7)] + 7AV • V 

- -V ■ [A(V k u x V k2 - V k u 2 V kl )} 

= -AV[trVu l2 (u) - (V kUl V k2 - V k u 2 V kl )} + 7AV • V 

- 1aV[(W + 2) 7 )]. 



Thus, by (12.91) . we can easily deduce that 

i V ■ [(VA + AV) 
AVj + 7AV • V = 7AV • V 



--V ■ [(VA + AV)9] + -V ■ (VuV - VVu T ) (3.26) 
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On the other hand, it follows from (l2"77j) . O and fl^TU]) that 

2 A ( VA + AV W - ^A{VuV - VVu T ) + jV (3.27) 
= -~/AA{I + V) - 7/ - ^*tV[-cj 12 (u) 

+ 7 ]J - -^VmF - VVu T ) 
= -jAA(I + V) - ~[-trVw 12 {u) + j{trV + 2)]I 

- -A(VuV - VVu T ) 
= - 7 AA(I + V)- -(VaUiV^ - V k u 2 V kl )I 

+ -A(v«y + y Vu T ) 

= -jAA(I + V). 

Hence we conclude from fl3.25l) - fl3.27l) that 

d t [V -V -A(I + V)V6] + u ■ V[V • V - A(I + V)V9] (3.28) 
= [u(u) + -yA][V-V- A(I + V)V9\. 

which concludes the proof of the lemma, since the above quantity will remain zero all 
the time with zero initial data. 

Remark 3.4. In order to carry out the mechanical background of the above lemma, 
we again go back to the deformation tensor F. By its definition, the following identity 

V • F T = 

is proved in [TH [T5] . Insert (12.11) into the above equality, we find 

diRij + RkjdiVik + VikdiRkj = 

which implies 

diVu + RijdiRij + RijVikdiRkj = 

since R is an orthogonal matrix. But then employing (12. 2p a straightforward calculation 
yields ( Kim . 



4 Proof of Theorem 2.2 



In this section we prove Theorem 12.21 Weak dissipations on the strain V and the 
gradient of rotation angle V# are found by introducing auxiliary functions w and 
below. The way of defining such functions reveals the intrinsic dissipative nature of 
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the system. In what follows, || • || is used to denote the standard L 2 (R 2 ) norm, < •, • > 
and (■, ■) the R d and L 2 (R 2 ) d inner products with d — 1, 2 or 4. The proof is divided 
into four steps. 

Step 1. Standard Energy Estimates. 

By taking the L 2 inner product of the second and third equations of system ( 11.2ft 
with u and V, and then adding up the resulting equations together, we obtain 



f (\u\ 2 + \V\ 2 )dx + ti\\Vu\\ 2 
2dtJ R 2 

f „\u\ 2 + \V\ 2 
= - u- V— — — 1 — — h <u,Vp> dx 
Jb 2 2 

+ f < V, ~(VuV + Wu T ) + -u 12 {u)(VA - AV) > dx 
Jb? 2 2 

+ / < m,2V • V > + < V, D(u)> dx 
Jb 2 



B? 



+ / < u, V ■ (VV )> dx- <V, -j(VA - AV) > dx. 
Jb 2 Jb 2 2 

On the other hand, the fourth equation of system ( II. 2ft gives 

-— [ 2trV dx = - [ u-VtrV dx + [ trD(u) dx 

2 dt J R 2 J R 2 J R 2 

+ / tr{yuV + \{VA - AV)(ui 2 {u) - 7)) dx. 
Jb 2 2 

Adding up the above two equation together, we have 

li [ (M 2 + \V\ 2 + 2trV) dx + fi\\Vu\\ 2 (4.1) 
2 dt J R 2 

f \ u \2,\y\2 r 

= - tt-V( +trV)+ < u, Vp > dx+ trD{u) dx 

Jb 2 2 J R 2 

+ / < V,\{VuV + VVu T ) + \u l 2{u){VA- AV) > dx 
Jb 2 2 2 

+ / < u, 2V • V > + < V, D(u) > dx + / tr(VuV) dx 
Jb 2 Jb 2 

+ / < u, V • (VV T ) > dx - [ <V, \i{VA - AV) > dx 
Jb 2 Jb 2 2 

+ f t r ±((VA-AV)(u Vi (u)->y)) dx. 
Jb 2 z 

Noting that V is symmetric, and with the aid of the first equation of ( 11.21) . it is 
rather easy to see that 

- J R2 u ■ \{\u\ 2 + \V\ 2 + 2trV)+ <u,Vp> dx + J R2 trD{u) dx = 0, 
J R2 < u, 2V • V > + < V, D(u) > dx + J R2 tr(VuV) dx = 0. 
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Noting the definition in (12.61) and (12. 9p . one has 

/ <V, \{^uV + Wu T ) + \u l2 {u){VA - AV) > dx (4.3) 

+ [ < u, V • ( W T ) > dx- [ <V, \i{VA - AV) > 

+ / t4((VA-Ay)(o;i 2 (u)-7)) dar 
Jr? 2 

= / < V, VuV > + < u, V ■ ( W T ) > dx 
Jr? 

+ / -(wi 2 (u) --f)(tr(VA- AV)+ < V,VA-AV > ) dx 
Jr 2 2 

= 0. 

Then, combining (14. ip . (14. 2 j) with (14.31) . we have 

-- / (|w| 2 + | V| 2 + 2trV) dx + Vwf = 0, 

Finally, by using (13.21) . we arrive at 

n^, / (M 2 + l^u -^ 2 | 2 + |^i2 + ^ 21 | 2 ) dx + /i||Vu|| 2 = 0. (4.4) 
2 at ,/ R 2 

By taking the L 2 inner product of the second and third equations of system (II. 2p 
with u and 2V, and then adding up the resulting equations together, we obtain 



Ydt 



2 

u\ l + 2\V\ l dx + I u-V{— + \V\ 2 )dx + {u,Vp) 
r? Jr 2 2 



= (V, VuV + VVu T + wi 2 (w)(VM - AV) - j(VA - AV) 
+ (u, V • (VV^)) + (m, 2V • V) + (2V, D(u)) + /i(u, Au), 

which gives that 

\j t J R H 2 + 2 \ V \ 2 dx + ^\\Wu\\ 2 (4.5) 

< I|V«||l~||V|| 2 

< C||W||||AV||(||AVu|| + ||Vu||) 

< C/i||Vy|| L 2(||Vn|| 2 + ||AV«|| 2 + ^||Ay||). 

/i 

In order to exploit the higher order energy estimates, we apply A to the second 
and third equations of system (ll.2p . and then take L 2 inner product of the resulting 
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equations with Au and 2AV, respectively, to yield 



1 d 
2di 



Au\ 2 + 2\AV\ 2 dx + (Au,VAp) (4.6) 



i? 2 

2 



+ / u-V{^^ + \AV\ 2 )dx 

R 2 2 



(AV, A [VuV + VVu T + u 12 (u) (VA - AV) 

- i(VA - AV)]) + (Au, AV ■ (VV T )) + (Au, 2AV ■ V) 

+ (2AV, AD(u)) + fi(Au, AAu) - (Au, [A(u ■ Vu) - u ■ VAw] 

- (2AV, [A(u ■ W) - u ■ VAV] ) . 



Applying similar arguments as before, we conclude that 

f R2 u ■ Vi}^f + \AV\ 2 )dx + (An, VAp) 

- fi(Au, AAu) = fi\\AVu\\ 2 , (4.7) 
(Am, 2AV • V) + (2AV, AD(u)) = 0. 

To estimate the rest terms of the right side of (j4.6j) . and also for later use, we 
need the following Lemmas. Their proofs can found in many literatures, see [H], for 
example. 

Lemma 4.1. Assume v G W k,2 (R 2 ), k > 3. The following interpolation inequalities 
hold. 



1. Forl<s< k, 



\\v\W < C|H| 1-5J ||V a v 
||Vw|| L 4 < CllwH 1-2 ^!) || V s Vf 1| 2T^FTT, 

||Av|| L 4 < C||t)|| 1_ ^2)||V s Av||^2). 



2. For2<s< k, 



\\v\\ L oc < C||v|| 1_ = ||V s n||-, 

< CH^II 1 a + 1 1| V s Vi>|| s+1 . 

The following Lemma will also be used for several times below. 

Lemma 4.2. Assume f,gE H S (R 2 ) for some s > 0. Then there exists a universal 
constant C depending only on s such that 

l|v s (/(?)|| <c(||/|| L oo||v^|| + |M|l~||V7II). 



16 



Now we estimate the second line of (14. 6p . Invoking the definition of 7 in (12.91) . and 

(4.8) 



with the aid of Lemma I4.2[ we deduce that 

(AV, A [VuV + VVu T + u 12 (u)(VA - AV) - j(VA - AV)] 
< C\\AV\\(\\AV\\\\Vu\\ L *o + ||VAw||||y||ioo + II VAu|| \\V\\ 2 Loo 



+ II Vnll^o 



l« II Ay || + ||y| 



1 



IA- 



1 



2 + trV 



< C f i\\V\\ H 2(\\Vu\\ 2 + ||VAu|| 2 , „ 

/i 2 



IAy| 



provided that 

which will be verified below. 

Similarly, it is easy to show that 



< 1, 



Au, AV ■ (VV T )) I 
< II Vllxoo II AVII II VAxtl 



(4.9) 



(4.10) 



< C/i||y||^ 2 (||VA W || 2 + — ||Ay| 



At last, it remains to estimate the last line of (14.61) . Again, by Sobolev imbedding 
theorem and Lemma 14.21 it leads to 

- (Au, [A(u -Vu)-u- VAu]) - (2AV, [A(u ■ VV) - u ■ VAV^j (4.11) 
< C||A«||(||V«||ioo||A«|| + ||VVtt||||tt||z,°=) 

+ c||Ay||(||v«||x»||Ay|| + ||vvu||||y|| r «.) 



< C{1 + ii){\\u\\ m + ||y||^)(||V«|| 2 + ||VA«|| 2 + ^||Ay|| 2 ). 

Combining all these estimates (I4.7p - (l4.1ip with (14.51) . and then noting (I4.6p . we can 
finally conclude that 



1 d 



u 



R 2 



< C{l + n){\\u 



+ |Aw| 2 + 2|Vf + 2| AV\ 2 )dx + /i(|| Vw|| 2 + || VAu|| 2 ) (4.12) 

1 



H 2 



y 



H 2 



VAu 



\AV\\ 2 ) 



provided (14.91) holds. 

Step 2. Weak Dissipation for The Strain Matrix V. 

To get dissipative energy estimates, by (I4.12p . it is clear that we should explore the 
dissipation of V. 
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Let 



w = Au + -V • V. 



(4.13) 



Now we apply A to the second equation of system i\1.2L and respectively, -V- to 
the third one, and then adding up the outcomes to get a new equation. In terms of w 
introduced in f)4.13p . this new equation can be written as 



where 



w t + u ■ Vw + VAp = /iAu> H — Au + /. 



(4.14) 



/ = - [A(u -Vu)-u- VAu] - - [V • (u ■ VV) 

A 4 



(4.15) 



a 



V(V • V)] + AV • (VV T ) + -V • [{VuV + VV 

/i 



u 



+ u) 12 (u)(VA - AV) - i{VA - AV)] . 

Take the inner product of f)4.14p with w in (L 2 (R 2 )) 2 , and use the similar arguments 
as in (14.21) . we have 



1 d 



\w\ 2 dx + n\\ Vw\\ 2 = (Ap, V ■ w) + -(Au, w) + (f, w 



(4.16) 



As before, we should estimate the right side of (14.161) term by term. Firstly, it is 
obvious that 

-|(Ait,io)| < tIIV^II 2 + \\\Vu\\ 2 . (4.17) 
/i 4 /j, 6 

On the other hand, noting that V ■ u = 0, we can use integration by parts and 
Lemma 14.21 to estimate the last term of the right side of (I4.16P as follows: 

\(f,w)\ < |([A(^-) -u t Au,],V t w J >)\ + \(A(VV T ),Vw)\ (4.18) 



1 

+ - 

Vro 



(j(VwV + Wu T ) + lu 12 (u)(VA - AV) - ry(VA - AV)] . 



+ - 



[Vjt ■ (uiV jk ) - UiiVkVjk)] , ViW, 



< C\\Vw\ 



Itill^oo || Au\ 



1 



I Vw|| L oo || Vw| 
1 



+ -11^11^11^11+ \\V \\ L °o AV 

<c(p + -)(\\u\\ H * + \\v\M 
[I 

x (||Vw|| 2 + ||AVm|| 2 + || Vwf + -j ||AV|| 2 ) 

/i 



V|||oc||Vu| 



18 



provided ||V||ioo < 1. 

It also remains to estimate the first term of the right side of (14.161) . For this purpose, 
we apply the divergence operator for the momentum equation of system ( 11.21) to get 

Ap = -tr(VuVu) + V ■ [V ■ (VV T )} + 2V ■ (V • V). (4.19) 

Thus, by ( I4.19p . Lemma H~Tl Lemma H~2l and Sobolev imbedding theorem, we can easily 
deduce that 

| (Ap, Vic) | < | (tr ( Vu Vu) , Vic) | (4.20) 
+ (V [V- (VV T )}, Vw) +2|(V- (V-V),Vw)\ 

<C\\Vw\\(\\Vu\\% + \\V\\ L ~\\AV\ 

+ ^||V W || 2 + ^||V-(V-V)|| 2 
< C(l + fi)(\\u\\ H 2 + \\V\\ H 2)(\\Vu\\ 2 + \\VAu\\ 2 

+ iiv^f + -4l|Ay|| 2 ) + ^||VH| 2 + -liv ■ (v ■ v)\\ 2 . 

M II ^211 II / 4H II ^11 V ^11 

We must deal with the last term of the above estimates. Now Lemma 13.31 enters the 
argument in this stage. Noting the definition of V x , it is easy to see that V ■ V ± = 0. 
Thus, by Lemma 14.11 and a straightforward calculation, we have 

-||V ■ (V ■ V)\\ 2 = -\\V • (AVV9)\\ 2 (4.21) 
fi u 

< -(l|v^||^l|v^|| L4 ) 2 + (llFlUooiiv^n) 2 

<-||V^||||V^||||Ay||||A^|| + ^H^lli. ||A^|| 2 

U U 



< C(l + [i){\V\\ m + |MM(||VAu|| 2 + — ||Afl|| 2 + — ||AV|| 2 ) 



u 2 u 

provided (14.91) is satisfied and 



||V0|| < 1 (4.22) 

holds. 

Finally, the combination of (14. 161) - (14. 2 II) leads to 

— / \w\ 2 dx + u\\Vw\\ 2 (4.23) 
dtJ R 2 

< C{u+-){\\u\\ H 2 + \\V\\ H 2)(\\Vu\\ 2 + ||AVm|| 2 
+ II VH| 2 + \\\^V\\ 2 + \\\A9\\ 2 ) + -JVu\\ 2 
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provided (|Q21) holds. 

Step 3. Weak Dissipation for V9. 

At this stage, to extract the dissipative nature of the system, it is clear that we 
need also explore some kind of weak dissipation of Ad. For this purpose, we introduce 
the auxiliary variable G: 

= Au+-V ± 9. (4.24) 

Rewrite the momentum equation as follows: 

u t + u ■ Vu + Vp = fj.Au + 2V ± 9 + V ■ (VV T ) + 2AVV6. (4.25) 

Then, apply A to (14.251) and ^V- to the last equation of (11. 2ft . and then add up the 
resulting equations together, we obtain a new equation in terms of O: 

Q t + u ■ V6 + VAp = A6 + - Aw + g. (4.26) 

A* 

where 

g = AV • (VV T ) + 2A(AVV6) + -V ± 7 (4.27) 

/i 



[A(«.V«)-«.VA«]-l(-^^) 



Similarly, we take the L 2 inner product of (14.261) with to get 

1 4- / |e| 2 rfx + ^||vef = (Ap,v-e) + -(Au,e) + (p,e). (4.28) 

2 at J R 2 /i 

By the expression of in (14.241) and the first equation of system (jl.2p . it is easy to see 
that 

(Ap,V-6) = 0. (4.29) 
A straightforward computation implies 



-|(Au,e)| < ^||ve|| 2 + ^||v«|| 2 . (4.30) 

fi 1 1 4 fj, 6 



To estimate the last term in (14.281) . we insert (14.271) into (g, 0) and then use inte- 
gration by parts to yield 

\(g, 6)| < \([A(u iUj ) - Ui Au 3 ] , VA-) | + ^ ( (^f^f ) , ©) | (4-31) 

+ |(A(w T ),ve)| + -|( 7 ,v ± -e)| + 2|(A(Ayve),e)|. 
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By a similar argument as before, we have 

[A(u iUj ) - UiAuj] , ViBj) + | (A(VV T ) } V6) | + 2| (A(AW0), 6) | (4.32) 



< C||VG|| [(||u|| L oc||A«|| + ||Vu|Ucc||V«||) + ||F|U«||AF|| 

+ c||ve||(||y||^||A^|| + ||vy|| L4 ||v^|U4) 

< c(i + ^)(\\u\\ H , + \\v\\ H 2 + -\\ve\\) 

n 

x (||Vu|| 2 + HAV^II 2 + ||V0|| 2 + -JAVf + - 2 \\A6\\ 2 ), 

/i /i 

where we used Lemma I4TT1 and Lemma |4~2"1 On the other hand, by f)2.9p . it follows that 

2 



Vi^Vim J 



(T.v-L-e)! 



(4.33) 



< ^(IIVTiiiiiveii^iieii^ + nviix-llVtiiiveii) 

< ^(||Vri||||Vfl||i||e||3||Afl||i||ve||^ 

< C(l + i)(||e|| + ||F|| H2 + i||V0||) 

/i /i 

x (||V^|| 2 +||V0|| 2 + ^||A^|| 2 ). 

A 4 

The combination of (I4.3ip - fl4.33p implies that 



||V«||V0||) 



1 



\(g,e)\<C(^ + -)(\\u\\ H2 
/i 



1 



fla + He|| + -||ve||) 

Vu\\ 2 + ||AV M || 2 + ||V6|| 2 + -A\A6\\ 2 + -AAV\\ 2 ). 

II* j2 z 



(4.34) 



Finally, combining (l^T28l) - fl^T30l) with fT4~3"4"l) . and noting the definition of 9 in ffl~24"l) . 

we arrive at 



d_ 

It 

< ^\\vu\\ 2 + c(i, 2 + -)(\\a\\ ll , -i \\v\\ ir ^ 1 1 c-> 1 1 ) 



\e\ 2 dx + ve|| 2 

B? 

C „_ „ 2 ; Is 

/'" A*' 



x 



(||V«|| 2 +||AV M || 2 + ||V0|| 2 + ^||Ay|| 2 ). 



Step 4. A Priori Weakly Dissipative Energy Estimates. 
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This step ties everything together. Firstly, by Hodge's decomposition, we calculate 

AV = VV-V-VxVxV (4.35) 

fv 2 (VxV 12 -v 2 v n ) -Vx(v 1 v 12 -v 2 v u )\ 
\y 2 (v 1 v 22 -v 2 v 21 ) -Vi(ViV 22 - v 2 y 21 y 

= VV • 1/ + (V^trV - V X (AV • V). 

Thus, noting (13.21) . and by (14. 13ft and Lemma I4~2l we obtain 

\\AV\\ < Cn(\\AVu\\ + \\Vw\\ + ||AtrV||) (4.36) 
< Cn(\\AVu\\ + IIVHI) + C^||Ay||||y||roo 

which implies that 

||AV|| < C»(\\AVu\\ + IIVHI) (4.37) 

provided that 

L oo < (4.38) 



(4.39) 



2fiC 

On the other hand, by (14.241) . 



U\ve\\<c^\\Au\\ + 

\||A0||<C7/i(||AVu|| + ||Ve||). 



Now we can proceed with the dissipative energy method. Plug (I4.39j) into (14.121) . 
we find that 

~l {\u\ 2 + \Au\ 2 + 2\V\ 2 + 2\AV\ 2 )dx + ^{\\Vu\\ 2 +\\VAu\\ 2 ) (4.40) 
2 at J R 2 

< C(l + /i)(||u|| H2 + \\V\\ H .)(\\\7u\\ 2 + HVA^II 2 + ||VH| 2 ). 

Similarly, insert (14 .3 7B and (14.391) into (14. 23ft and (14.351) . and then add the resulting 
inequalities together, we obtain that 

4/ (kr+|e| 2 )dx + /i(||V^|| 2 ||+V0|| 2 ) (4.41) 
< C{f + i)(||«||^ + UVllfl* + ||0||)(|| W|| 2 + || AV^f 



+ ||V«;|| 2 +||Ve|| 2 )+ — ||V^|| 2 . 



c 
(1 

Combining (I4.40H with (14.411) . we finally arrive at 
d 



\(\u\ 2 + \Au\ 2 + 2\V\ 2 + 2\AV\ 2 ) + \w\ 2 + \Q\ 2 }dx (4.42) 
at J R 2 L J 



+ /i(|| vu\\ 2 + \\vAuf + 1| vhi 2 + II ve|| 2 ) 

^||Vn|| 2 + C( M 2 + - 

M d /I 
x (||vu|| 2 + ||VAu|| 2 + ||VH| 2 + liven 2 ). 



< -iivh 2 + c{f + -)(\\u\\ H , + \\v\\ H * + 1|6||) 
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Therefore, if the initial data is sufficiently small, we can find some positive time T, 
such that 

11*11*° + ^ + » e » £ 2C(fb) <443) 
for all < t < T. Moreover, in this case, (14.421) implies that 



d 
It 



[ \(\u\ 2 + \Au\ 2 + 2\V\ 2 + 2\AV\ 2 ) + \w\ 2 + | O | 



dx 



^(||v«|| 2 + ||vam|| 2 + ||vhi 2 + livef) < — i|Vw| 



Thus, noting (I4.13p . (I4.24p . and by integrating the above inequality from to t, < 
t <T, with respect to t, we obtain 

(ll«ll^ + ll^ll^ + lkl| 2 + H0|| 2 ) (4-44) 

POO 

+ m / (IIVm||h2 + ||Vwf + ||ve|| 2 )^ 

Jo 

< c(ll«oll^ + IIHII^ + IKII 2 + l|e || 2 ) + - / ||v«|| 2 dt 

A 4 Jo 

< C(l + ^(Huoll^ + \\V \\ 2 H2 + ||V0 o ||hi), 

/i 

where in the last inequality, we used the basic energy law (14.41) . Then, it follows from 
(I4.44p that if the initial data satisfies 

\\u \\ 2 H2 + \\V \\ 2 H2 + \\V9\\ 2 m < M{ f + (4.45) 

for a big enough constant M (M > 2C 3 ) independent of /i, then (14.431) will still hold 
for all the latter time t > T with a uniform constant C independent of t and /i. 
Furthermore, when (14.451) holds, we can deduce from (14.441) that 

(NI^ + IMI^ + IMI 2 + ll©l| 2 ) (4-46) 

/■oo 

+ A* / (l|Vu|||a + \\Vw\\ 2 + ||V6|| 2 )dA 
Jo 

< C(\\u \\ 2 H2 + \\V \\ 2 H2 + IKf + ||e || 2 ) + — / ||Vu|| 2 dt 

A* Jo 

< ^(1 + Zi 4 ) 



C 2 (l+/i 10 ) 

provided C is big enough. 

At last, we should show (l4~9L (@~22]) and fjOSfl . By flOSl and Sobolev imbedding 
theorem, we have 

c/j, 2 

\\V\\l°° < c M Ur 2 < 



2C(/i 3 + 1) 
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for appropriate constants c and C. Thus, we have 



. 1 < 1, for a > 1, 
1 1 < i, for /i < 1, 



which means that (14.91) and (14.381) hold. To show (14.221) . we deduce from (I4.43P and 

IIV.II < 2MIIBH + HW < < 1. 

The proof of Theorem 2.2 is completed. 
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